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INTRODUCTION 


Several  recent  studies  have  found  a  lack  of  coherence  between  ground 
motion  waveforms  recorded  at  stations  separated  by  only  a  few  tens  of  meters. 
Reference  1  covers  the  station-to-station  waveform  coherence  for  near-field 
explosion  accelerograms  recorded  on  a  ni.ie-station  array  at  Pahute  Mesa, 

Nevada  Test  Site.  For  this  array,  at  a  range  of  6  km  from  a  5.6  Mj_  under¬ 
ground  explosion  with  an  interstation  spacing  of  100  meters,  strong  incoherent 
signals  were  found  above  5  Hz  on  all  components.  This  incoherence  was  attrib¬ 
uted  to  scattering.  Reference  2  discusses  earthquake  seismogram  coherence 
for  a  nine-station  array  with  an  interstation  spacing  of  CO  m.  This  tempo¬ 
rary  array  was  located  near  the  Pi  non  Flat  Observatory  in  California. 

Several  events  with  magnitudes  between  3.0  and  4.2  and  hypocentral  distances 
ranging  from  10  to  50  km  were  analyzed.  For  this  data  set,  P  waves  were 
found  to  be  coherent  up  to  the  25  to  30-Hz  region,  and  S  waves  up  to  about 
15  Hz.  This  variability  in  ground  motion  over  such  short  distances  raises 
questions  about  the  appropriateness  of  using  single  station  point  measure¬ 
ments  of  the  ground  motion  field  for  the  inversion  of  source  functions  as 
well  as  the  use  of  the  peak  ground  acceleration  (PGA)  statistic  based  on 
single-point  measurements. 

Reference  3  analyses  the  near-field  motions  from  the  1979  Imperial 
Valley  Earthquake.  This  data  set  was  obtained  from  a  closely  spaced  array 
of  five  stations.  This  study  confirmed  that  incoherence  of  higher  frequency 
ground  motions  could,  because  of  averaging  effects,  reduce  the  high-frequency 
input  to  buildings  with  foundations  of  large  areal  extent.  Analysis  of  this 
particular  data  set  suggested  that  base  averaging  reduction  factors  of  10 
to  30  percent  would  have  occurred  for  foundations  50  m  in  diameter  for  fre¬ 
quencies  above  20  Hz.  No  significant  reduction  would  have  occurred  for 
frequencies  below  5  Hz. 

This  frequency-dependent  spatial  variability  in  ground  motion  has  also 
been  observed  on  some  small-scale,  high-explosive  tests  conducted  in  dry 
alluvium  (Ref.  4).  An  array  of  six  triaxial  accelerometers  spaced  at  six 
azimuths  at  the  20-m  range  from  a  buried  5-lb*  detonation  yielded  waveforms 

*To  convert  from  lb  to  kg,  multi-'y  by  4.535924  E-01. 
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which  were  incoherent  above  30  to  35  Hz.  In  an  effort  to  determine  whether 
these  variations  were  induced  by  instrumentation,  source  asymmetry  or  scat¬ 
tering  due  to  geologic  variability,  a  series  of  small-scale,  high-explosive 
experiments  were  performed.  This  paper  describes  the  results  of  these 
experiments,  an  effort  to  develop  a  simple  seismic  survey  technique  which 
could  be  used  to  determine  the  coherence  function  of  the  geological  structure 
at  a  given  test  site,  and  the  physical  interpretation  of  the  data. 
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THE  INITIAL  EXPERIMENT— THE  ART  2  TEST 

The  Array  Test  Series  (ARTS)  was  a  series  of  single  and  multiple  5-lb 
explosive  tests  conducted  in  dry  alluvium.  The  ARTS  was  formulated,  designed 
and  implemented  to  investigate  the  quantitiative  effects  of  finite  explosive 
sources  upon  observed  motion  fields.  The  single-burst  tests  in  the  series 
were  intended  to  thoroughly  characterize  the  single-burst  ground  motion 
environment  (Ref.  5).  One  of  these  tests,  ART  2  (Table  1),  was  intended  to 
determine  the  degree  of  azimuthal  variation  in  ground  motion  levels  resulting 
from  a  single  charge  event.  The  experimental  layout  for  this  test  is  shown 
in  Fig.  1.  Six  triaxial  servo  accelerometers  were  placed  at  the  20-m  range 
at  six  different  azimuths.  These  accelerometers  were  recorded  by  3-channel 
digital  event  recorders  at  a  rate  of  200  samples/second  with  a  5-pole  low-pass 
Butterworth  filter  at  70  Hz.  The  original  accelerograms,  along  with  low-pass 
(30-Hz  corner)  and  high-pass  (40-Hz  corner)  filtered  versions,  are  shown  in 
Figs.  2  and  3.  Fourier  acceleration  spectra  for  the  waveforms  are  shown  in 
Fig.  4.  Examination  of  the  acceleration  traces  and  the  spectra  in  these 
figures  reveals  considerable  similarity  in  the  various  waveforms  below  30  Hz, 
while  large  azimuthal  amplitude  variations  are  observed  for  frequencies  above 
30  to  35  Hz.  The  accelerograms  show  that  the  high-frequency  energy  arrives 
after  the  initial  coherent  low-frequency  response,  with  a  great  variability 
in  wave  shape  as  well  as  amplitude.  The  vertical  and  radial  spectral  shapes 
(Fig.  4)  are  similar  between  5  and  35  Hz;  above  35  Hz,  variations  up  to 
15-20  dB  occur. 


3 


I 


AFWL-TR-86-133 


TABLE  1.  Description  of  experiments. 


Experiment 
ART  1 


ART  19 


ART  II 


ART  1 1 1-2 


Shotgun  Tests 


Description 

Single  5-lb  charge  buried  at  1  m.  Recording  array  -  6 
tri axial  accelerometers  at  different  azimuths  at  the 
20-m  range. 

"Huddle  Test."  Single  5-lb  charge  buried  at  1  m 
recorded  by  6  triaxial  accelerometers  in  a  closely 
spaced  group  at  the  20-m  range. 

The  "Pit  Snots."  Series  of  five  5-lb  buried  charges 
all  fired  at  the  same  ground  zero  and  recorded  by  6 
triaxial  accelerometers  spaced  at  equal  azimuths  at 
the  20-m  range.  The  initial  charge  was  detonated  in 
situ;  subsequent  charges  were  fired  in  a  sand-filled 
pit. 

Single  5-lb  buried  charge  recorded  by  6  accelerometers 
at  the  10-m  range  and  6  accelerometers  at  the  30-m 
range.  Accelerometers  spaced  at  equal  azimuths. 

Seisgun  recorded  by  12  or  24  vertical  geophones  at 
ranges  from  10  to  30  m. 
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RAINED  SAND 


ACCELEROMETER 
3  COMPONENT 


5  lbs  CHARGE 

Experimental  layout  for  the  ART  2  Test  Event. 
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Figure  2.  Filtered  vertical  acceleration  records  from  ART  2. 
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Figure  3.  Filtered  radial  acceleration  records  from  ART  2. 
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EXPERIMENTAL  DESIGN 

Three  possible  sources  for  the  observed  azimuthal  variations  in  wave¬ 
forms  have  been  identified:  (1)  nonuniformity  in  instrumentation  response, 
(2)  asymmetry  in  the  explosive  source,  and  (3)  geologic  inhomogeneity. 

In  order  to  experimentally  investigate  the  first  possibility-instrumental 
irregularity— the  ART  19  test  (Table  1),  also  called  the  "huddle"  test,  was 
performed.  In  this  test,  all  seven  3-component  accelerometers  were  placed 
in  a  closely  spaced  group  at  the  20-m  range  (Fig.  5)  from  a  5-lb  charge 
detonated  at  a  depth  of  1  m.  Individual  vertical  and  radial  spectra  from 
this  test  are  compared  in  Fig.  6.  The  resulting  spectral  scatter,  as 
observed  in  these  figures,  is  at  most  3  to  4  dB,  over  the  band  from  5  to 
70  Hz— much  less  than  that  observed  for  the  azimuthal  data.  In  fact,  the 
higher  frequency  scatter  (40  to  70-Hz  region)  shown  in  these  figures  may  not 
reflect  true  instrumental  irregularity,  since,  as  seen  by  examination  of 
Figs.  5  and  6,  the  accelerometers  which  are  physically  closest  together  in 
the  field  do  possess  the  greatest  degree  of  similarity  in  spectral  shape. 

An  additional  set  of  experiments,  the  ART  II  series  (Table  1),  was 
designed  to  determine  whether  azimuthal  irregularities  in  ground  motion 
result  primarily  from  propagation  path  variations  and  scattering  by  geologic 
inhomogeneity  or  explosive  source  asymmetry.  The  test-bed  layout  for  the 
ART  II  tests  is  shown  in  Fig.  7.  Ground  motions  were  recorded  by  tri axial 
servo  accelerometers  and  3-channel  digital  event  recorders  at  six  equally 
spaced  azimuths  at  a  range  of  20  m.  The  ART  II  series  of  tests  consisted  of 
five  separate  5-lb  detonatlons—all  fired  at  the  same  test-bed  with  the  accel¬ 
erometers  remaining  in  the  same  position  from  shot  to  shot.  The  initial 
detonation  was  fired  at  a  depth  of  1  m  in  situ.  The  resulting  crater  was 
then  excavated  and  a  new  5-lb  charge  emplaced  in  the  resulting  pit  with 
uniform  sand  rained  around  the  charge.  This  process  was  repeated  four  times. 
Thus,  the  ART  II  series  of  tests  consisted  of  five  shots,  one  in  situ  and 
four  pit  shots—all  using  the  same  test-bed  and  an  identical  recording  array. 

The  goal  of  the  ART  II  tests  was  to  determine  whether  source  assymetry 
or  geologic  factors  were  the  dominant  cause  of  the  observed  scatter.  By 
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Figure  5.  Experimental  layout  for  Huddle  Test. 
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Figure  7.  Test-bed  layout  for  the  ART  II  tests. 
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repeating  the  experiment  five  times  at  the  same  location,  five  sets  of  tri- 
axial  accelerograms  and  spectra  were  obtained  for  each  azimuth.  If  the 
observed  data  scatter  were  produced  by  geologic  inhomogeneity,  then  one 
would  expect  very  similar  acceleration  waveforms,  as  well  as  amplitude 
spectra,  to  be  observed  from  shot  to  shot  at  the  same  recording  station. 

The  detonation  of  a  5-lb  charge  is,  of  course,  not  strictly  a  repeatable 
phenomenon,  since  the  initial  in  sit!»  detonation  destroys  a  portion  of  the 
test-bed  by  formation  of  the  crater.  Excavation  of  the  severely  disturbed 
material  comprising  the  crater  and  subsequent  emplacement  of  the  charge 
within  the  pit  was  intended  to  make  the  experiment  as  repeatable  as  possible. 
The  transition  from  the  in  situ  test  to  the  pit  sequence  was  the  most  severe. 
The  in  situ  shot  produced  a  crater  roughly  2.2  m  in  diameter.  The  crater  was 
excavated  producing  a  pit  about  3  m  in  diameter  and  1.3  m  deep.  The  dimen¬ 
sions  of  both  crater  and  pit  increased  slightly  with  each  shot.  The  final 
pit  was  about  3.5  m  in  diameter  and  the  final  crater  about  the  same  diameter. 

While  the  ART  2  and  the  ART  Phase  II  (the  pit  shots)  experiments  investi¬ 
gated  azimuthal  ground  motion  variations  at  the  20-m  range,  the  ART  1 1 1 - 1 
experiment  (Tab.  1)  was  designed  to  investigate  the  influence  of  range  upon 
the  observed  azimuthal  variations.  If  scattering  by  geologic  inhomogeneity 
is  the  dominant  factor  producing  the  azimuthal  variation,  then  one  might 
expect  the  degree  of  variation  to  increase  with  range  of  the  receiver  array. 
The  ART  III-l  test  consisted  of  a  5-lb  cylindrical  charge  detonated  in  situ 
at  a  depth  of  1  m.  The  ground  motions  were  recorded  at  12  locations  using 
triaxial  servo  accelerometers  and  3-channel  digital  event  recorders.  Six 
recording  stations  were  placed  at  equal  azimuths  at  the  10-m  range,  and  the 
remaining  six  were  placed  at  equal  azimuths  at  the  30-m  range. 


13 


jin*.  *u.  ■ 


z 


HA  Wl  JW^  ATWJWMl ^  W*  l  W!Ai\njl  jui  f^f  ■  #V»l  A-.1  i*./i  Hi  , ..^»  .  **.-.  mu  1  **.-  »..-.  mr-  ^.-i  an«.-.  arll 

AFWL-TR-86-133 

STOCHASTIC/DETERMINISTIC  CLASSIFICATION  TOOLS 
AND  DATA  ANALYSIS 

One  approach  to  quantify  the  observed  variation  in  the  data  sets  is  to 
compute  the  coherency  between  station  pairs  in  a  manner  analogous  to  that 
employed  in  Ref.  1.  The  coherency  is  defined  by 

y  x  y 

where  <j>  xy(co)  is  the  cross  power  spectrum  and  $  xx(oj)  •  $  yy(w)  the  power 
spectra  of  the  two  time  series. 

Smoothing  must  be  applied  to  the  SDectra  prior  to  estimation  of  the 
coherency  in  order  to  minimize  variances  at  each  frequency.  For  the  coherency 
estimates  computed  here,  a  4-point  lag  window  was  applied  to  the  spectra  prior 
to  estimation  of  the  coherency  factors. 

Coherencies  were  computed  between  the  0°  data  and  each  of  the  other 
recording  stations.  The  time  windows  used  were  128  samples,  or  0.64  s  in 
length.  This  length  window  encompasses  all  of  the  wave  train.  At  this  rela¬ 
tively  short  range,  it  is  effectively  impossible  to  separate  the  various  com¬ 
ponents  of  the  waveform. 

The  computed  coherencies  for  the  station  pairs  are  shown  in  Fig.  8. 

The  coherence  factors  for  the  vertical  and  radial  components  are  higher  than 
might  be  intuitively  expected  after  examination  of  the  spectral  plots  in 
Fig.  4.  However  as  noted  in  Ref.  3,  the  coherency  will  be  unity  when  the 
signals  at  two  stations  are  related  by  any  linear  transfer  function,  so  the 
coherence  function  is  not  a  direct  estimator  of  the  variation  of  amplitude 
across  an  array. 

Based  upon  the  results  of  the  coherency  analysis,  an  alternate  approach 
for  analysis  of  the  data  was  taken.  The  set  of  amplitude  spectra  of  the 
azimuthal  observations  was  treated  as  a  statistical  ensemble.  Figure  9  gives 

14 


AFWL-TR-86-133 


COHERENCE  WITH  O*  DATA 


VERTICAL  TRANSVERSE  RADIAL 


Frequency 


Figure  8.  Computed  coherencies  for  ART  2. 


an  example  of  this  technique  applied  to  the  ART  2  data.  The  mean  and  variance 
of  the  spectral  ensemble  were  calculated  as  shown,  assuming  a  lognormal  dis¬ 
tribution.  These  figures  illustrate  y\  the  mean,  and  y  ±  o,  the  mean  ±1 
standard  deviation  for  the  vertical,  radial,  and  transverse  data  from  ART  2. 

The  vertical  and  radial  data  show  small  standard  deviations  between  5  and  30  Hz. 
Beyond  35  Hz,  the  width  of  the  y  ±  o  deviation  band  increases  until,  at  40  Hz, 
the  expected  scatter  in  the  Fourier  modulus  is  8  dB  or  greater.  The  trans¬ 
verse  motions,  on  the  other  hand,  show  at  least  8-dB  scatter  throughout  the 
entire  frequency  band;  although  the  overall  mean  amplitude  is  reduced  about 
10  dB  from  the  radial  and  vertical  amplitude  levels.  The  lack  of  coherence 
below  5  Hz  on  all  components  has  several  possible  sources.  This  effect  was 
first  thought  to  be  due  to  signal -induced  tilts  in  the  accelerometers.  Inte¬ 
gration  of  the  acceleration  traces  and  subsequent  attempts  to  correct  the 
resulting  velocity  and  displacement  records  did  not  significantly  change  the 
estimates.  Aliasing  and/or  noise  may  be  the  cause  of  this  low-freouency 
incoherence.  The  amplitude  level  in  this  region  is  very  low  relative  to  the 
higher  frequencies.  Even  a  small  amount  of  aliased  energy  would  significantly 
affect  the  apparent  coherence  in  this  region  which  is  well  below  the  corner 
frequency  of  the  5-lb  source.  At  higher  frequencies  the  aliased  energy  level 
is  insignificant  relative  to  the  amplitude  of  the  direct  source  energy. 

Another  measure  of  the  data  scatter  in  the  frequency  domain  is  the 
coefficient  of  variation  (CV)  (Ref.  6).  The  CV  is  the  ratio  of  the  standard 
deviation  to  the  mean,  CV  =  a/y.  The  CVs  for  the  ART  2  spectra  are  given  in 
Fig.  10.  For  ART  2,  the  CV  on  the  radial  and  vertical  components  increases 
with  frequency  (CV  <  0.5  from  5  to  30  Hz);  the  transverse  component  CV  remains 
approximately  constant  at  a  relatively  high  level  (>  0.5)  across  all  frequen¬ 
cies,  suggesting  that  the  transverse  component  is  uncorrelated  at  all 
frequencies.  By  way  of  contrast,  1 ig.  11  shows  the  very  low  values  of  the  CV 
between  5  and  60  Hz  for  the  radial  and  vertical  components  cf  the  Huddle 
(ART  19)  test.  The  advantage  of  the  CV  as  an  estimator  of  the  data  scatter 
across  the  array  is  that  all  stations  in  the  array  may  be  used  in  the  single 
estimate,  as  opposed  to  only  a  pair  of  stations  in  the  coherence  estimate. 

Amplitude  spectra  of  the  vertical  acceleration  waveforms  recorded  at  the 
0°  azimuth  for  all  five  of  the  ART  II  pit  shots  are  compared  in  Fig.  12. 
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Figure  10.  Spectral  coefficient  of 
variation  curves  for  the 
ART  II  pit  shots  compared 
with  ART  2. 
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Figure  11.  Spectral  coefficient  of  variation  curves  for  the  vertical 
and  radial  components  of  the  Huddle  Test. 
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Figure  12.  Comparison  of  vertical  acceleration  spectra  for  the  five 
shots  of  the  ART  II  series. 
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The  sets  of  spectra  recorded  at  the  other  azimuths  and  the  other  components 
are  similar  in  character  to  the  0°  vertical  set.  The  four  pit  shot  spectra 
are  quite  similar  in  shape  and  amplitude,  while  the  in  situ  shot  spectrum 
exhibits  a  shape  which  is  similar  to  that  of  the  pit  shot  spectra,  but  higher 
in  amplitude  for  almost  all  frequencies.  Two  of  the  pit  shots  are  almost 
identical  in  the  frequency  domain;  pit  shots  3  and  4  also  group  closely 
together  over  most  of  the  frequency  band.  The  greatest  change  in  amplitude 
from  shot  to  shot  occurs  in  the  30  to  35-Hz  region.  This  area  of  the  spectrum 
is  fairly  flat  for  the  in  situ  shot;  however,  a  spectral  hole  develops  for  the 
pit  shots  and  deepens  with  each  successive  pit  shot.  Even  though  the  shape 
of  the  source  spectrum  may  change  from  shot  to  shot  (which  is  to  be  expected 
since  introduction  of  the  backfilled  cavity  results  in  more  energy  lost  in 
cratering,  greater  attenuation  and  a  larger  cavity),  the  CV  should  change  very 
little  if  the  source  remains  cylindrically  symmetric  for  all  shots.  Figure  10 
shows  the  azimuthal  coefficient  of  variation  plotted  as  a  function  of  frequency 
for  the  vertical,  radial  and  transverse  components,  respectively,  for  all  five 
shots.  The  normalized  standard  deviation  curves  are  very  similar  for  all  five 
shots  out  to  the  system  anti-alias  filter  at  70  Hz.  The  coefficient  of  varia¬ 
tion  curves  from  ART  2,  which  was  conducted  near  the  pit  shot  site  but  on  a 
different  test-bed,  are  also  plotted  in  this  figure.  The  change  in  the 
coefficient  of  variation  curves  between  the  two  sites  is  much  larger  than  the 
variation  from  shot  to  shot  for  the  pit  shot  sequence.  The  relatively  small 
scatter  in  the  coefficient  of  variation  curves  for  the  five  events  of  the  pit 
shot  sequence  (the  in  situ  shot  plus  the  four  pit  shots)  suggests  that  source 
asymmetry  plays  an  unimportant  role  in  the  development  of  the  azimuthal  varia¬ 
tion  in  ground  motion  response. 

The  ART  1 1 1-2  (Table  1)  shot  was  designed  to  determine  the  degree  to 
which  azimuthal  variation  increases  with  range  from  the  source.  In  this  test, 
ground  motions  were  recorded  by  circular  arrays  at  both  the  10-m  and  30-m 
ranges.  Coefficient  of  variation  plots  for  the  vertical,  radial  and  transverse 
spectra  for  both  the  10-m  and  30-m  ranges  are  given  in  Fig.  13,  These  results 
support  the  contention  that  the  coefficient  of  variation  increases  with  range. 
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NEW  SITE  CHARACTERIZATION  PROCEDURE 

One  goal  of  this  study  is  the  development  of  a  standardized  survey  tech¬ 
nique  which  can  characterize  the  degree  of  geologic  inhomogeneity  present  and 
its  influence  upon  azimuthal  ground  motion  variations  across  a  site.  Such  a 
technique  requires  a  fairly  energetic  seismic  source  which  is  both  cylindri- 
cally  symmetric  and  repeatable.  The  Betsy  Seisgun* **  was  chosen  for  testing  as 
a  potentially  useful  source  for  these  sorts  of  surveys.  The  Betsy  device  is 
an  C-gauge  industrial  shotgun  converted  for  use  as  a  seismic  source.  The 
Betsy  fires  an  3-oz  slug  into  the  ground  surface,  providing  a  fairly  energetic 
seismic  source.  In  order  to  quantify  this  source,  a  series  of  tests  similar 
to  those  conducted  with  the  5-lb  charges  was  performed  with  the  Seisgun. 

Seisgun  tests  were  performed  within  the  ART  alluvial  test  site  at  two 
different  locations  about  500  m  apart.  One  of  these  sites  was  the  same  test¬ 
bed  which  was  subsequently  used  for  the  ART  III  10-m  and  30-m  array  tests. 

This  site  will  be  referred  to  as  the  ART  III  site;  the  other  site  will  be 
referred  to  as  Site  A.  All  Seisgun  tests  used  vertical  4.5-Hz  natural  fre¬ 
quency  geophones  as  motion  sensors.  The  geophones  were  recorded  by  a  24- 
channel  digital  seismic  recording  system  at  a  sample  rate  of  0.25  ms,  with  an 
anti -alias  filter  corner  at  1000  Hz. 

Two  separate  arrays  were  employed  at  Site  A.  One  array  was  composed  of 
12  geophones  spaced  at  30-deg  angles  at  the  20-m  range  from  the  Betsy.  The 
other  array  consisted  of  a  "Huddle"  test  in  which  all  12  geophones  were  placed 
in  a  closely  spaced  group  in  order  to  determine  the  uniformity  of  instrument 
response. 

One  objective  of  the  Seisgun  tests  at  Site  A  was  to  determine  the  repeat¬ 
ability  of  the  signal  from  shot  to  shot.  Figure  14  illustrates  the  mean  and 
variance  of  the  amplitude  spectrum  as  recorded  at  a  single  station  for  four 
successive  shots  of  the  gun  at  Site  A.  The  signal  is  quite  repeatable  from 

*Betsy  Seisgun  is  a  trademark  of  Betsy  Seisgun,  Inc,  Tulsa,  OK. 

**To  convert  from  oz  to  kg,  multiply  by  2.834  952  E-02. 
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Figure  14.  Mean  and  variance  of  spectra  for  four  shots  of  the  Selsgun 
at  the  same  site. 
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shot  to  shot.  These  shots  were  the  sixth,  seventh,  eighth,  and  ninth  firings 
into  the  same  hole  (i.e.,  the  gun  was  fired  five  times  prior  to  the  first 
recording  in  order  to  establish  proper  gain  levels  on  the  recording  system). 

The  signal  may  not  be  as  repeatable  for  the  first  few  shots  into  the  same  spot 
because  of  changing  source  characteristics  as  the  hole  dug  by  the  slug  becomes 
deeper. 

Figure  15  compares  the  coefficient  of  variation  of  the  amplitude  spectrum 
for  the  "Huddle"  test  with  the  coefficient  of  variation  computed  for  the  full 
12-channel  circular  array  test  at  Site  A.  The  full  array  coefficient  of 
variation  is  perhaps  0.3  to  0.4  greater  above  35  Hz  than  that  observed  in  the 
"Huddle"  test,  again  implying  that  the  observed  azimuthal  variations  in  ampli¬ 
tude  levels  for  the  full  circle  array  test  at  Site  A  are  not  due  to  instrumental 
irregularity. 

A  number  of  tests  were  conducted  with  the  Seisgun  at  the  ART  III  site. 
Recordings  were  made  at  the  10-m,  20-m,  and  30-m  ranges.  Twelve  vertical 
geophones  were  positioned  at  30-deg  increments  at  each  range.  Recordings  were 
performed  for  several  shots  of  the  gun  with  a  total  of  24  geophones  at  the 
10-m  and  20-m  ranges.  The  10-m  phones  were  then  shifted  to  the  30-m  range  and 
several  more  shots  were  recorded.  Objectives  of  these  gun  tests  were:  (1)  to 
determine  repeatability  of  array  coefficient  of  variation  data  for  several 
shots  of  the  Seisgun,  both  into  the  same  hole  and  when  the  position  of  the 
Seisqun  is  shifted  slightly  from  the  original  shotpoint,  and  (2)  to  determine 
how  the  coefficient  of  variation  changes  with  increasing  range  (unfortunately, 
we  were  not  able  to  record  unclipped  waveforms  at  the  10-m  range). 

Figure  16  contains  a  comparison  of  the  20-m  coefficient  of  variation  for 
three  Seisgun  shots  at  the  ART  III  site  and  one  shot  at  Site  A.  The  coeffi¬ 
cient  of  variation  curves  for  the  three  shots  at  the  ART  III  site  are  quite 
similar,  particularly  at  the  lower  frequencies.  The  Site  A  curve  diverges 
fairly  widely  from  the  ART  III  curve  above  30  to  40  Hz  (at  60  Hz  the  Site  A 
curve  approaches  a  CV  value  of  1.0  while  the  ART  III  CV  is  near  0.5),  indi¬ 
cating  the  difference  in  degree  of  geoloqic  scattering  between  the  two  sites. 
Two  of  the  ART  III  curves  in  Fig.  16  are  from  the  7th  and  22nd  shots  fired 
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Figure  16.  Comparison  of  the  coefficient  of  variation  curves  for  three 
Seisgun  shots  at  the  ART  1 1 1-2  site  and  one  shot  at  Site  A. 
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into  the  same  hole;  the  remaining  curve  was  obtained  from  the  4th  shot  fired 
into  a  point  on  the  surface  which  was  displaced  approximately  0.5  m  from  the 
other  shotpoint.  The  excellent  agreement  between  the  curves  obtained  from  the 
two  different  shotpoints  shows  that  the  observed  azimuthal  variations  are 
indeed  a  result  of  geologic  scattering  and  not  source  irregularity  or  asym¬ 
metry.  In  Fig.  17,  a  comparison  is  made  between  the  coefficient  of  variation 
curves  obtained  from  the  20-m  and  30-m  circle  arrays  at  the  ART  1 1 1-2  site. 

As  in  the  ART  III  5-lb  CV  comparison  (Fig.  13),  the  CV  values  for  the  larger 
range  circle  array  are  higher  over  most  of  the  frequency  band  than  the  close- 
in  array  values.  This  increase  of  azimuthal  variability  with  increasing  range 
from  the  source  is  what  would  be  expected  for  scattering  due  to  geologic 
inhomogeneity. 
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RELATION  OF  OBSERVED  AZIMUTHAL  VARIABILITY  TO  TEST  SITE  GEOLOGY 

The  tests  discussed  In  this  report  were  conducted  at  the  same  test  site 
In  the  Rio  Grande  Valley  south  of  Albuquerque,  New  Mexico.  The  subsurface 
geology  consists  of  dry  alluvium  down  to  the  water  table  at  a  depth  of  approxi¬ 
mately  75  m.  The  site  lies  in  a  small  shallow  basin  which  in  earlier  times 
contained  a  playa.  In  general,  the  site  is  underlain  by  unconsolidated  eolian 
sands,  alluvium  and  lacustrine  deposits  15  to  30  m  thick.  Beneath  this  mate¬ 
rial  are  the  tertiary  Santa  Fe  and  Galisteo  formations,  180  to  300  m  thick 
(Ref.  7).  A  typical  near-surface  P  wave  velocity  section  is  shown  in  Fig.  18 
(Ref.  8). 

Near-surface  deposits  are  fairly  variable  with  intermittent  caliche  beds 
present.  Trenching  at  other  test  sites  in  the  area  has  revealed  discontinuous 
subsurface  caliche  beds  typically  on  the  order  of  0.5  to  1.0  m  in  thickness 
(Ref.  8).  In  an  effort  to  quantify  the  nature  of  the  subsurface  variability 
at  the  ART  2  test  site,  18  boreholes  were  drilled  within  the  confines  of  the 
test-bed  (Fig.  19).  Each  borehole  was  drilled  to  a  depth  of  6.1  m.  Standard 
penetration  tests  in  each  hole  determine  the  number  of  hammer  blows  required 
to  drive  a  sampling  tube  a  unit  distance.  This  blow  count  is  related  to  the 
in  situ  density  and  compressive  strength  of  the  subsurface  material  (Ref.  9). 
Figure  20  displays  the  blow  counts  per  meter  for  each  borehole  dvjwn  to  the 
6-m  depth.  Figure  19  contains  isopach  maps  showing  the  lateral  distribution 
of  the  areas  of  similar  blow  counts  as  a  function  of  depth. 

A  thorough  understanding  of  the  physical  processes  responsible  for  the 
observed  azimuthal  variability  in  the  ART  2  test  requires  that  a  quantitative 
relationship  be  established  between  the  subsurface  material  property  informa¬ 
tion  and  the  observed  ground  motion  characteristics.  Aki  and  Richards 
(Ref.  10)  present  such  a  relationship.  Following  Chernov  (Ref.  11),  they 
characterize  the  statistical  nature  of  geologic  media  by  the  autocorrelation 
of  the  velocity  fluctuation. 

For  u  =  -6C/C0,  where  C0  =  velocity,  the  normalized  autocorrelation  function 
is  defined  by 
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boreholes  within  the  ART  2  test-bed  and  isopach  map  of  blow  counts. 
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Figure  20.  Plots  of  blow  counts  vs  depth  (data  points  from  all  18  bore 
holes  are  included). 
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N(r)  =  <u(rJ)  •  u(r'+r)>/<y2> 


where  r  denotes  oosition. 


Two  statistical  distributions  are  considered:  (1)  the  case  where  N(r)=  e 
(exponential  distribution),  and  (2)  the  Gaussian  distribution  given  by 
N(r)  =  e~r  ^a.  The  quantity  a_  is  the  measure  of  scale  length  of  the  medium 
inhcmogeneity  called  the  correlation  distance. 


- iri/a 


Once  a  medium  is  described  statistically  in  terms  of  N(r),  then  the  Born 
scattering  approximation  may  be  used  to  determine  the  strength  of  the  scatter¬ 
ing  as  a  function  of  frequency  for  the  medium.  Following  Ref.  10,  the  ratio 
of  scattered  energy  to  total  energy  is  given  as  follows: 

si/I  =  8<u2>  •  k4  •  a3  •  L  for  ka  <  1 
and 

51/ I  =  2<u2>  •  k2  •  a  •  L  for  ka  >  1. 


where  51/ I  is  the  ratio  of  scattered  energy  to  the  total  energy  carried  by  the 
wave,  <u2>  is  the  average  of  the  square  of  the  velocity  perturbation,  K  is  the 
wave  number,  a  is  the  scale  length  or  correlation  distance  of  the  medium,  and 
L  is  the  length  of  the  travel  path. 

Making  the  assumption  that  there  is  a  linear  correlation  between  the  blow 
count  value  at  a  particular  point  and  the  P  wave  velocity  at  the  same  point, 
we  shall  attempt  to  relate  the  observed  azimuthal  variability  in  the  ART  2 
data  with  the  blow  count  data  using  the  scattering  theory  discussed  above. 

The  mean  was  computed  for  the  set  of  blow  count  data  for  all  of  the  holes. 

This  yields  a  mean  value  of  approximately  50  blows/0.30  m.  This  value  was 
treated  as  C0  in  the  scattering  equations.  Therefore,  the  autocorrelation 
function  of  the  following  ratio  was  then  computed:  (BCi-BC)/BC,  where  BCi 
is  the  individual  blow  count  value,  and  BC  is  the  mean  value  of  all  blow 
counts.  The  autocorrelation  functions  were  computed  separately  for  each  hole. 
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Values  of  the  autocorrelation  functions  of  the  blow  count  data  for  all  of 
the  holes  are  plotted  in  Fig.  21.  Also  plotted  in  Fig.  21  are  several  theo¬ 
retical  autocorrelation  functions  for  different  values  of  the  correlation 
length  a,  for  both  the  Gaussian  and  exponential  distributions.  Based  on  exami¬ 
nation  of  these  plots,  the  theoretical  distribution  function  which  best  fits 
the  data  is  apparently  an  exponential  function  with  a  scale  length  between 
2.0  and  3.0  m.  This  scale  length  applies,  of  course,  only  to  the  vertical 
distribution  of  inhomogeneity.  The  horizontal  inhomogeneity  scale  length 
could  well  be  different  in  character.  The  horizontal  spacing  of  the  boreholes 
(  =  6.7  m)  is  much  too  wide  to  allow  any  sort  of  reasonable  estimate  of  the 
horizontal  scale  length  (this  spacing  yields  a  nyquist  wavelength  of  13.4  m) . 

Using  the  Born  equations,  the  ratio  of  scattered  energy  to  total  energy, 
(<5I/I),  was  computed  for  the  ART  2  test-bed  using  several  values  of  the  inhomo¬ 
geneity  scale  length  as  shown  in  the  plots  in  Fig.  22.  The  ka«l  case  is 
appropriate  for  most  of  the  frequency-scale  length  combinations  under  con¬ 
sideration  here.  The  ka»l  case  is  appropriate  for  scale  lengths  greater  than 
2  m  in  the  higher  frequency  region  (>20  Hz).  Both  curves  indicate  that,  for 
scale  lengths  of  2.0  to  3.0  m,  significant  (6 1/ 1  >  0.10)  scattering  begins 
near  10  Hz.  Significant  scattering  in  the  ART  2  test  data  occurs  above  30 
to  35  Hz  (Fig.  4).  There  are  several  possibilities  for  this  apparent  poor 
correlation:  (1)  the  Born  approximation  is  valid  only  for  small  values  of 
scattered  energy;  (2)  the  autocorrelation  function  which  was  derived  pertains 
only  to  the  vertical  distribution  of  inhomogeneity,  and  thus  the  scale  length 
which  is  appropriate  for  the  overall  propagation  medium  could  be  considerably 
different  from  2.0  to  3.0  m;  (3)  a  finer  spatial  sampling  of  the  near  surface 
blow  count  values  might  well  yield  a  smaller  scale  length;  and  (4)  the 
assumed  correlation  between  blow  count  data  and  P  wave  velocity  may  be 
incorrect. 


To  determine  the  effect  of  varying  the  scale  length  upon  the  theoretical 
scattered  energy,  6 1/ I  as  a  function  of  frequency  was  computed  for  correlation 
lengths  of  1.0  m,  0.5  n.  and  0.25  m  (Fig.  22).  For  these  scale  lengths,  the 
frequency  at  which  scattering  becomes  significant  is  in  the  25  to  30-Hz  region. 
This  result,  particularly  the  0.5  curve,  correlates  much  better  with  the 
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EXPONENTIAL  DISTRIBUTION 


GAUSSIAN  DISTRIBUTION 


Figure  21.  Experimental  autocorrelation  function  compared  with 
Gaussian  and  exponential  autocorrelation  functions. 

(The  lines  on  the  plot  are  the  theoretical  autocor¬ 
relation  functions;  the  experimental  autocorrelation 
functions  from  all  holes  are  represented  by  data  points.; 
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Figure  22.  Plot  of  scattered  energy  as  a  function  of  fre¬ 
quency  for  different  scale  lengths. 
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observed  threshold  of  incoherence  in  the  ART  2  test.  A  correlation  scale 
length  in  the  neighborhood  of  0.5  to  1.0  m  is  also  not  unreasonable  in  light 
of  the  results  of  the  "Huddle"  test  (ART  19)  in  which  all  of  the  instrumenta¬ 
tion  was  placed  in  a  group  approximately  1  m  wide  (Fig.  5).  Examination  of 
the  spectra  from  this  test  in  Fig.  6  reveals  that  the  spectra  of  the  endmost 
instruments  separated  by  a  distance  of  approximately  1  m  begin  to  diverge; 
whereas,  the  spectra  from  instruments  spaced  more  closely  correlate  very  well. 
Shallow  trenching  in  the  general  test  site  area  suggests  that  a  typical  caliche 
bed  thickness  is  on  the  order  of  0.5  to  1.0  m,  although  a  typical  lateral 
extent  may  be  slightly  greater  (Ref.  8). 

This  simple  forward  model  indicates  that  the  azimuthal  variability  in  the 
ART  2  data  can  be  explained  with  a  scattering  model  based  upon  the  Born 
approximation.  The  implication  is  that  the  particular  standard  penetration 
test  array  did  not  properly  characterize  the  appropriate  scale  length  due  to 
Inadequate  horizontal  or  vertical  spatial  sampling,  since  blow  count  values 
were  obtained  by  counting  the  number  of  blows  required  to  drive  the  sampling 
tube  a  distance  of  0.6  m.  This  sample  interval  yields  a  Nyquist  wavelength 
of  1.2  m. 
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DISCUSSION  AND  CONCLUSIONS 

The  results  of  the  ART  experiments  combined  with  the  Betsy  Seisgun 
observations  have  shown,  at  least  for  this  particular  alluvial  test  site,  that 
the  predominant  cause  of  azimuthal  variation  in  ground  motion  response  is 
inhomogeneity  in  near-surface  geologic  material  rather  than  very  near-source 
asymmetry.  Considering  that  this  particular  site  in  general  is  thought  to 
have  a  fairly  uniform  near-surface  geologic  composition,  the  degree  of  azimuthal 
variation  observed  from  these  small  tests  is  fairly  remarkable.  For  the  short 
ranges  involved  here,  the  8-gauge  shotgun  provides  a  repeatable  and  symmetric 
source  to  be  used  in  the  development  of  a  standardized  seismic  survey  technique 
for  characterization  of  these  azimuthal  variations.  Surveys  of  this  type 
should  be  performed  at  several  sites  to  further  constrain  the  relationship 
between  the  spectral  coefficient  of  variation  and  the  variation  in  the  geologic 
structure.  Such  surveys  will  aid  greatly  in  the  interpretation  and  under¬ 
standing  of  ground  motion  data  from  small-scale  field  tests.  For  a  small- 
scale  test  performed  at  the  ART  2  site,  deterministic  interpretations  of  the 
data  cannot  be  made  above  the  30  to  35-Hz  region.  Coefficient  of  variation 
surveys  could  also  place  limitations  upon  the  interpretation  of  near  source, 
regional,  and  teleseismic  measurements  of  earthquakes  and  explosions  made  at 
a  single  point.  These  surveys  may  be  useful  in  explaining  the  lack  of  coher¬ 
ence  between  individual  station  pairs  in  small  seismic  arrays,  such  as  the 
Pinon  Flat  array  (Ref.  2). 

In  this  study,  an  initial  attempt  was  made  at  relating  the  available 
subsurface  information  at  the  ART  2  site  with  the  observed  waveform  variation 
in  the  test.  This  attempt  was  not  very  successful,  probably  due  to  insuffi¬ 
cient  spatial  resolution  of  the  subsurface  sampling  performed.  Forward 
modeling  argues  that  the  scale  lengths  which  are  present  in  the  ART  2  test-bed 
are  perhaps  on  the  order  of  0.5  to  1.0  m  rather  than  2.0  to  3.0  m  as  deter¬ 
mined  from  the  autocorrelation  of  the  standard  penetration  test  blow  count 
data.  Use  of  other  sampling  techniques,  such  as  the  seismic  cone  penetrometer 
test  (in  which  the  force  necessary  to  drive  a  rod  into  the  ground  is  measured 
and  a  geophone  in  the  rod  tip  is  used  to  perform  a  velocity  survey),  could 
obtain  much  finer  spatial  resolution  both  horizontally  and  vertically, 
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allowing  a  greatly  improved  characterization  of  the  random  properties  of  the 
medium. 
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